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Abstract 

We study the possible distortion of phase-space in the decays 
K — > TTTTTT, which may result from final state interference among the 
decay products. Such distortion may influence the values of slope 
parameters extracted from the Dalitz plot distribution of these decays. 
We comment on the consequences on the magnitude of violation of the 
I AI 1= 1/2 rule in these decays. 
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1 Introduction 



In most cases, strangeness changing hadronic decays obey the | A/ |= 1/2 
rule. Experimental information on the cases when the rule is violated may 
help in understanding its origin. K-decays provide an important laboratory 
where one could learn more about this important issue In particular, 
by measuring the slope terms of the K ^ Sn Dalitz plot, as suggested by 
Weinberg one can probe the | A/ |= 1/2 rule with good precision. 

The slope parameters of the Dalitz plot distribution in non leptonic K^^^ 
decays were measured some time ago with limited precision More ac- 
curate experimental results have been published recently |^ attracting 
renewed theoretical interest 

The Dalitz plot distribution for — > Svr can be parametrized by a series 
expansion of the form [Q] 

\ M \'^ocl + gX + hX'^ + kY'^ + ... (1) 



with 



and 
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So — ^ Si — [Pk — Pi) 

where pi are the four momenta of the pions {i = 1,2,3) and the label 3 
denotes the odd pion in a decay. 

The coefficients g, h and k are not available from theory, but can be extracted 
from the experimental Dalitz distribution. In the literature there are other 
definitions of slope parameters. We choose a definition compatible with the 
Particle Data Group p|. 



2 Bose-Einstein interference 

Bose-Einstein Correlation (BEC) was observed for the first time by Gold- 
haber, Goldhaber, Lee and Pais |^ in pp interactions. They found that the 
angular distribution of like-sign pion pairs was different from that of unlike- 
sign pairs. In general terms, the effect can be understood as the tendency of 
identical bosons to occupy the same phase space. As a consequence, identical 
bosons are correlated in their momenta, opening angles distributions, etc. 
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The impact of BEC on the Dahtz plot distribution in hadronic decays of 
charm mesons has been studied before 0. It has also been demonstrated 
to play an important role in other phenomena M in high energy physics, 



and on the extraction of standard-model parameters [0, 10|. More recently, 



the effect of Bose-Einstein interference on measurement of collective flow has 



been considered 111 



In order for Bose-Einstein correlation to be present among identical par- 
ticles in the final state of any reaction, two aspects must be fulfilled: i) the 
source of pions must be finite, ii) the emission must take place chaotically to 
a certain degree. 

In the cases under study i.e., the decays: 



TT+TT^TT^ (2) 

K+ 7r-7r+7r+ (3) 

7r±7r°7r° (4) 

and 

K° ^ ttOttOttO (5) 

condition i) is satisfied by the need of form factors to describe the decay, and 
the well known finite size of the K mesons ||12|| . 

In hadronic decays, such as those under study, particles are produced par- 
tially through hadronization in a not a entirely coherent process. The inco- 
herence although small, we believe it to be sufficient to satisfy condition ii), 
requiered for the presence of correlation. Figure 1 illustrates the decay. The 
bubble represents the space-time region where particles are produced. 

Correlations among the decay products of a particle and the pions pro- 
duced in the main reaction, as described in Ref. §], will not be considered 
here. We will rather regard the decay itself as a particle production process 
in which interference may arise. 

The BEC are commonly described in terms of a two-particle correlation 
function: 
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where P{pi,p2) is the joint probabihty amphtude for the emission of two 
bosons with momenta pi and p2, and P{pi) and P{p2) are the single produc- 
tion probabihties. 



The Bose-Einstein correlation among identical mesons has been used to 
probe the space-time structure of the intermediate state right before hadrons 
appear |^ |13], Q in high energy elementary particle and nuclear collisions. 



One parametrizes the effect assuming a set of point-like sources emmiting 
bosons. These point like sources are distributed according to a density p{r). 
The correlation function can then be written as, 

R{Pi,P2) = J P(n)p(^"2) I i'BEiPi,P2) P d^ri£r2, (7) 

where pi,P2 are the momenta of the two bosons, i/jbe represents the Bose- 
Einstein symmetrized wave function of the boson system, with Jy p{r)d^r = 
1. 

Taking plane waves to describe the bosons, one obtains the correlation 
function for an incoherent source: 

i?(pi,P2) = l+|^(p(r-)) |2, (8) 

where J-'{p) represents the Fourier transform of the density function p{r). 

Phenomenological parametrizations of the effect have been proposed, to 
describe the quantum interference during the hadronization in high energy 



reactions. For a recent review see Ref. JT^. Here, we use the GGLP 
parametrization to gain an idea of the impact of particle interference on the 
slope parameter. In a more detailed study [[15[, we will address other possible 
parametrizations. 

The GGLP parametrization is one of the most commonly used. It is given 

by: 

i?(p-l,P2) = l + Ae-^«', (9) 

where is the Lorentz invariant, = —{pi — ^2)^, which can be written 
also as, = ~ 4m^; here mi2 is the invariant mass of the two pions 
and m^r the mass of the pion under consideration. The parameter A lies 
between and 1, and reflects the degree of coherence in the production. The 
radius of the source is defined by r = hc^/J3[f'm]. The presence of Bose- 
Einstein correlations will modify not only the invariant mass spectrum of 
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like-charged but also that of unlike- charged pions. This reflection of BEC, 
known as residual correlation, has been studied so as to make sure that the 
reference sample of unlike-charged pions used to subtract the effect from the 
like-charged pions is free from any other correlations. 

In particle production processes with much phase space high multiplicity, 
residual correlations tend to be minimal. Nevertheless, some studies claim 
that this is not always the case, and that using unlike-charged pions as a 
reference may be questionable. 

In a particle production process yielding where only three particles {e.g., 
K — > tttttt) residual correlations are expected to be significant. 

3 Simulation of Bose-Einstein correlations 

In this letter we want to estimate the effects of BEC on the phase space of a 
three body decay. In particular, on the phase space of the decay K — > tttttt. 
We will take the approach used in Ref. 0, where the BEC are simulated 
by simply weighting each event. The MC generator consists of a three-body 
decay, with appropiate masses for the decay products. Each event is weighted 
according to: 

W = X[{l + Xe~^'^^, (10) 

where the product is taken over all pairs («, j) of like-charged pions. 

For cases when only two like-charged pions are present in the final state, as 

in the reactions (2-4), the product reduces to 

W = l + Xe-^'^\ (11) 
However, for the decay (5) it becomes 

W = {l + Ae-^^i2)(i + Xe-f^'^^il + Ae"'^^"). (12) 

where the Qij of all possible pairs are taken into account. We performed sim- 
ulations for different values of (3 and A. The coherence parameter A controls 
the strength of the effect, but does not modify the shape of the distortion. 
The decay is not necessarily a completely chaotic process, and A does not 
have to be 1. In fact, if hadronization did not take place during the de- 
cay, one would expect a completely coherent process in which A = 0, and 
Bose-Einstein interference would not be present. Hadronization in the decay 
introduces some degree of incoherence giving rise to A values between and 
1. The exact value would be obtained by fitting the correlation function, as 
in the case of the source radius. On the other hand, there may be production 
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mechanisms other than hadronization that introduce some degree of incoher- 
ence. The study of BEC in particular decays may help to disentangle the 
production processes involved. 



4 Effects on the Dalitz plot distribution 

Figure 2 shows the Dalitz plot distribution and its projections on X and 
Y, as defined in Eq. (1), for the decay 7r+7r'^7r°. The histograms 

show the projections with (full line) and without (crosses) BEC. The scat- 
ter plot is the ratio of the distribution with and without BEC, Where we 
took (3 = 25 GeV~^, which corresponds to r = 1 fm. Increasing (3 (Eq. 
(10)) would make the effect more marked. As mentioned above, the physical 
meaning of /3 = 4, 9, 16, 25 GeV~^ can be viewed in terms of the source ra- 
dius of r 0.4, 0.6, 0.8, 1.0 fm, respectively. The true value, however, should 
be extracted from experimental data measuring BEC in these kind of decays. 

In what follows, we will show the plots with (3 = 7 GeV~^, i.e., r ^ 0.5 
fm, which is motivated by the electromagnetic radius of the K~^. 

Ignoring electromagnetic interactions in the final state, the distributions 
for the modes (2), (3) and (4) would be identical. The correction for this ef- 
fect could be carried out dividing the phase space by the Gamow coefficients, 



as was, done in Ref. |T6|. Here, we try to focus on the effect of two-particle 
interference, and will not introduce this correction. In Ref. [|l^], we will 
present a complete anlysis, including the interplay of BEC and final-state 
electromagnetic interactions. 

Figure 3-5 show the amplitudes \ M vs. X, as obtained using PDG 
values (Table 1) for the different decay processes. The band represents 
the uncertainty in the parameters. The lines are the amplitudes one should 
observe, after correction for BEC, for = 7 GeV~^ and A = 0.25, 0.5, 0.75. 

The bands in Figures 3-5 show the errors according to Table 1 values. 
These curves are corrected and the corresponding errors are propagated. 
The corrected curves with errors are then fitted with the function of Eq. 1. 
In this way we mimic the procedure during data analysis and observe the 
change in the values of the parameters as they would be observed by the 
experimentalist. 

As one can see, the interference of identical particles in the final state changes 
the values of g and h. 

In the decay Kl tt^tt^tt^, all the pion pairs are weigthed according to 
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Decay 


mode 


g 


h 


k 


K+ ^ 




-0.2154±0.0035 


0.012±0.008 


-0.0101±0.0034 


K- 




-0.217±0.007 


0.010±0.006 


-0.0084±0.0019 


K^^ 




0.652±0.031 


0.057±0.018 


0.0197±0.0045±0.0029 


Kl^ 




0.678±0.008 


0.076±0.006 


0.0099±0.0015 


Kl^ 






-0.0033±0.0011±0.0007 
(-0.0061±0.0009±0.0005)t 





Table 1: Average values of Dalitz slope parameters for K ^ ^iti decays taken 
from the PDG 0. The NA48 collaboration 0] has recentely published the 
quadratic slope parameter for the decay K^j^ Tr'^Tr^vr^. The result is included 
here indicated with ^. 



Eq. (12). To good extension, the effect cancels out, and the Dalitz plot re- 
mains almost unchanged. Figure 6 shows the the amplitude using the PDG 
values (Table 1) for this decay. The scale in this plot has been expanded 
to make the effect visible. In this case the fit is done with a quadratic form 
without linear terms. 



5 Discussion and Conclusions 

The measured amplitud | M p as a function of X changes its shape once BEG 
among identical particles in the final state products is taken into account. 
The different curves shown in Fig. 4 represent the modified amplitude as a 
function of X for a pion source-radius of r = 0.5 fm, and a coherence factor 
of A = 0.25,0.5,0.75. The true values for (3 and A must be extracted from 
data. It may be that A and/or (3 are smaller than this, thereby reducing the 
impact on the amplitude. 

The Particle Data Group gives the average values for h and k (see 
Eq. (1)). These are shown in Table 1. 

After correcting the g and h parameters for BEG, by fitting 1 + gX + hX'^ 
to the corrected X distribution - shown in Figs. 3,4 and 5, one obtains the 
values given in Table 2. 

According to the | A/ |= 1/2 rule, one expects 

g{K^ ^ 7r^7r°7r°) ^ gjKj ^ 7r+7r-7r°) ^ _^ 



7 



Decay 


mode 


g (BEC corrected) 


h (BEC corrected) 


K+ ^ 




-0.1744±0.0027 


0.001±0.003 






-0.176±0.003 


-0.0008±0.0031 






0.692±0.005 


0.080±0.006 








-0.0052±0.0023 



Table 2: Central values of the BEC-corrected Dalitz slope parameters for 
^ 377 decays. The error is the result of the fit to the corrected curves 
once the errors are propagated in the correction procedure. 

and 

g{Kl ^ tt+tt-ttO) _ 

g{K+ 'k+'rH^) ■ ^ ' 

Using the PDG values of Table |l|, for g one obtains 

g{K^ 7r±7r°7rO) 



g[K — ^ TT TT TT 



^ tt+tt-ttOI 



n n n 



g{K+ tt+ttOttO) 
Using the values of Table these become 
g{K^ tt^ttOttO) 



3.0046 ±0.1726 



-3.1244 ±0.1073 



1.0399 ±0.0509 



9{K 
9{Kl 



n n n 



- tt+tt-ttO) 



9{K 

9{K_ 



n n n 



^ tt+tt-ttO) 



g{K^ 



-3.968 ±0.069 



-3.887 ±0.075 



0.980 ±0.014 



The ratios change by taking into account a phase space distortion due to 
the interference of identical particles in the final state. We have given only 
a general trend based on parameter values motivated by the electromagnetic 
structure of the kaon. A more precise statement would be possible once the 
experimentally extracted values of (3 and A are used in the analysis. One 
can see that the ratios given by Eq. (13) are very sensitive to a distortion of 
phase space, while that of Eq. (14) is not. Nevertheless although the size of 
the effect may be small, is clear that the observed violation of the | A J |= 1/2 
rule can be affected by BEC, and this should be taken into account in the 
anlaysis of data. 
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Figure Captions 



Fig. 1: Diagram representing the Kaon decay process. The bubble represents 
the space time region where pions are produced. 

Fig. 2: Dahtz plot distribution and its projections for the decay tt^tt^tt^. 
The histograms represent the distribution with (full line) and without 
(crosses) a contribution from BEC. The distributions have been nor- 
malized to I 1= 1 at X = and y = 0. The scatter plot shows the 
ratio of the distribution with and without BEC with parameters of the 
BEC simulation being (3 = 25 GeV~^ and A = 0.5. 

Fig. 3: Amplitude | M p vs. X for the decay tt+tt+tt", as obtained 

from experiment (shaded band) and the central value of the am- 
plitude that one would obtain (solid, shaded and dotted, lines) after 
correction for BEC effects with A = 0.25, 0.5, 0.75. The band represents 
the experimental uncertainties on the average values. 

Fig. 4: Amplitude | M p vs. X for the decay K~ — > 7r^7r^7r+, as obtained 
from experiment (shaded band) and the central value of the am- 
plitude that one would obtain (solid, shaded and dotted, lines) after 
correction for BEC effects with A = 0.25, 0.5, 0.75. The band represents 
the experimental uncertainties on the average values. 

Fig. 5: Amplitude | M p vs. X for the decay -k^-k^-k^ ., as obtained from 

experiment (shaded band) and the central value of the amplitude one 
would obtain (solid, shaded and dotted, lines) after correction for BEC 
effects with A = 0.25, 0.5, 0.75. The band represents the experimental 
uncertainties on the average values. 

Fig. 6: Amplitude | M p vs. X for the decay tt^tt^it^, as obtained from 

experiment (shaded band) and the central value of the amplitude one 
would obtain (solid, shaded and dotted, lines) after correction for BEC 
effects with A = 0.25,0.5,0.75. The band represents the experimental 
uncertainties on the average values. 
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